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Second-order Supersymmetric Operators and Excited States 

Micheal S. BergeiO and Nail S. Ussembaye\|j 
Department of Physics, Indiana University, Bloomington, IN 4^405, USA 

Factorization of quantum mechanical Hamiltonians has been a useful technique for some time. 
This procedure has been given an elegant description by supersymmetric quantum mechanics, and 
the subject has become well-developed. We demonstrate that the existence of raising and lowering 
operators for the harmonic oscillator (and many other potentials) can be extended to their super- 
symmetric partners. The use of double supersymmetry (or a factorization chain) is used to obtain 
non-singular isospectral potentials, and the explicit expressions for the ladder operators, wave func- 
tions and probability densities are provided. This application avoids the technical complexities of 
the most general approaches, and requires relatively modest methods from supersymmetric quantum 
mechanics. 

PACS numbers: 03.65.Ge, 03.65.Fd, 03.65.Ca 

In supersymmetric quantum mechanics (SUSY QM) the form of the Schrodinger equation of a particle is determined 
by its Hamiltonian (h = 2m = 1): 

0\\ H = -d 2 x + V(x) 

(N . 

where V(x) is nonsingular on — oo < x < oo and has at least one bound state. If its ground state wave function tp 
and energy Eq are known, tl 
Oh linear differential operators: 



o 

(N 



r^ . and energy Eq are known, then it is always possible to factor the Hamiltonian H_ = H — Eq into a product of two 



a ' H m = -d 2 x + V W (x) = (-d x + W Q )(d x + Wo)= a\a 

^ | where Wo (x) is the superpotential related to the potential V_ (x) by the Riccati equation: 

Vf } (z) = W^x) ~ W^x) = $,'(aOM>(*). 

^ • The superpotential can be written in terms of the ground state wave function as Wq (x) = — ip' Q {x)/i()o{x) with primes 
^O , denoting the derivative. By reversing the order of A Q and Ao, one obtains a new Hamiltonian 

tff = A Q Al = -dl + Vi 0) (x) = -dl + Wi{x) + W&x) 

whose eigenvalues and eigenfunctions are related to those of H_ . For instance, the eigenfunctions of H + (H_ ) can 



o 
o 



be found by applying the operator Aq (Aq) to the eigenfunctions of H_ (H^_ '). Moreover, it turns out that these 

partner Hamiltonians have identical spectra (except for the ground state of H\ ) , so that having an exactly soluble 

problem for H_ guarantees derivation of the exact solution of H + [l| . 

K> , The generalization of this formalism to the excited states ip n (x) (n = 1,2, • • • ) of H_ is straightforward Q. It is 

j^ ■ done by introducing the operators 

Al = -d x + W n (x) and A n = d x + W n (x) 

with superpotentials W n (x) = — ij)' n {x) / ij) n {x) . Now we have two sets of potentials V± (x) = W%(x) ± W^{x) 
corresponding to the wave function %j) n (x). From the Schrodinger equation it follows that 



C3 



v[ n \x)=V^\x)+E Q -E n 

igenvalue corresj 
other hand, the potentials 



(where E„ is the energy eigenvalue corresponding to ip n {x)), so the potentials V_ (x) offer nothing new. On the 



Vl n \x)=2Wl{x)-V^\x) 
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are less trivial and more interesting. A few remarks about this generalization are in order here. Since 

H^[A n ^ k ] = A n A\A n ^ k = A n [HtUk] = (E k - E n )[A n ^ k ] 

one may conclude that A n ip k is an eigenfunction of the Hamiltonian H\_ with energy eigenvalue E k — E n (note that 
the normalization condition requires k > n). However, from the orthogonality of eigenfunctions ipn(x) we know that 
an excited state wave function changes sign at one or more points in space and consequently the superpotentials W n (x) 
are always singular for n > 0. Therefore, the proper behavior of the wave function V'fc(x) does not necessarily mean 
that A n ip k (x) = (d x + W n (x))ip k (x) is an acceptable solution, and the states of H_ are not, in general, in one-to-one 

(n) 

correspondence with the states of H\_ unless n = 0, which corresponds to the standard unbroken SUSY. That is, 
for excited states the degeneracy of energy levels typically breaks down and this breakdown explains relatively low 
interest in extensions of SUSY QM to the excited states of a potential. As an illustration consider the following 
exercise. One dimensional harmonic oscillator eigenfunctions are given by ip n (y) = H n (y)e~ v ' 2 where H n {y) are 
Hcrmitc polynomials and y = ^w/2x. The superpotentials are easily derived by using the well known recursion 
relations for Hermite polynomials: 

W n (y)= ^( H ^ {y) 



2 V H n (y) 
They can be used to generate the first two pairs of partner potentials according to V±{x) = W 2 (x) ± W^(x): 



* 0) (*) ~ - | and Vi%) = <*f + £; U«(x) = ^ -% and V?\x) = <*f - £ + J, 



r(0). 



The n = case is the standard isospectral SUSY QM (the lowest state of V_ (x) has no counterpart in the spectra 

+ 



of VJ. (x)), whereas the n = 1 case corresponds to the generalization of SUSY QM to the first excited state and one 



expects the breakdown of the degeneracy of energy levels. Indeed, the potential V+ (x) defined on the interval [0, oo) 
is known as isotonic oscillator [3[ . Its eigenfunctions 



Aiip„(y) = J- 



-r + y-- 
dy y 



e- y2/2 H n (y) 



are only nonsingular for odd n and consequently the energy level spacing is twice larger than that of its partner 
V_ (x), i.e. the energy levels are partially degenerate. In the next pair of potentials V_ (x) and V7 (x) degeneracy 
is completely absent |J|. Note that we investigate the degeneracy between the potentials V± for some fixed n and 

not between two unrelated pairs of potentials V± and V± . 

It is clear that in the traditional factorization method (i.e. the first-order SUSY QM) there is no way to use the 
excited states of the initial potential and at the same time avoid creating singularities in the partner potential. The 
second-order SUSY QM can be used to overcome this difficulty [5| . It can be formulated by the following intertwining 
relationship 

HA = AH 

where two different Hamiltonians H = —d 2 + V{x) and H — d 2 x + V(x) are intertwined by the operator of the second 
order in derivatives: 

A = dl+j](x)d x +j(x). 

There exist two kind of operators A - reducible and irreducible. If it is not possible to write A as a product of two 
first order operators with real superpotentials, then they are referred to as irreducible [6(. In this paper we generate 
families of strictly isospectral potentials by using reducible operators which are equivalent to two successive first-order 
intertwining relations with real superpotentials. In other words, the relationship between partner potentials in our 
framework is the same as that obtained by the confluent algorithm [7| , but the novelty of our approach is based on 
exploiting the non- uniqueness of factorization in the construction of partner Hamiltonians [8[ . In this case there is 
an atypical Hamiltonian _ff_j_" (singular for n > 0) that admits non-unique factorization. This Hamiltonian plays 
an intermediate role in the construction and we do not discuss the degeneracy of its energy levels - we merely use 
it to construct two isospectral potentials. Thus, we succeed in avoiding some common obstacles (such as singular 



potentials, failure of the degeneracy theorem that arises when dealing with excited states etc.) by using non-uniqueness 
of factorization at the intermediate stage. Continuing our exercise with the harmonic oscillator we show below how 
to construct two one-parameter class of potentials which are not the harmonic oscillator, yet have the same energy 
spectrum. We explicitly derive the associated wave functions and plot the probability densities. 
The starting point is to define the operators 

B n =8 X + f{x) + W n {x) and B\ = -d x + f(x) + W n {x) 
where f(x) is temporarily undetermined function. As usual we demand that HY 1 = BnB^ which results in 

H { + ] = -d 2 x + Vi n \x) + /'(.*) + 2W n (x)f(x) + f 2 (x). 
In order to arrive at the non-unique factorization H+ = H+ one should set 

f'(x)+2W n (x)f(x) + f 2 (x)=0 

and solve for f(x). This differential equation is known as the Bernoulli equation (a specific example of the Riccati 
equation) and it possesses the following general solution which depends on an arbitrary integration constant C that 
can be considered as a free parameter: 



C + j*Jl{s)ds 

where C and xq are constants. A little calculation shows that the partner Hamiltonian H_ = B\B n is given by 

H { 1 1) = -d 2 x + Vi n) - 2f' n {x) = -dl + Vi n) (x). 

Obviously, the potential V+ {x) is free of singularities for n = and it is legitimate to discuss the degeneracy 
of energy levels of the potentials V_ {x), V_f_ ' (x) and V_ \x) (the latter is assumed to be nonsingular for certain 
values of the parameter C). However, as it was pointed out earlier for the excited states the potential V}_ n (x) has 
singularities corresponding to the zeros of the wave function ip n (x) . So while the potential V_ (x) is defined on the 

full real axis, the range of V+ (x) is composed of n + 1 disjoint regions and their energy spectra are non-degenerate 
or only partially degenerate as explained above. An interesting question arises naturally: is the degeneracy theorem 
still valid for the potentials V_ (x) and V_ (x) under the assumption that V_ (x) is free of singularities for certain 
C? We have found that the intermediate (or atypical) Hamiltonian Hi 1 admits two different factorizations due to its 
non-uniqueness, namely, AnA^ = B n B^. Hence, we can exploit the operators A n and B n together with their adjoint 
to answer the above question. 

The Schrodingcr equation H_ ipk = (Ek — E n )ipk implies 

H ( l l) [Bl A n 4, k ] = Bl L B n Bl L A n ip k = BlA n A\A n ^ k = (E k - E n )[BlA n ^ k ] 

where k ^ n, i.e. if ip k (x) is an cigenfunction of the Hamiltonian H_ with energy eigenvalue E k — E n , then B^A n ip k 

is an eigenfunction of H_ with the same energy. Similarly, from the Schrodinger equation H_ ip k = Ej^ 1 ip k 

(where in E[ , k denotes the energy level and (n) refers to the n th eigenfunction of the Hamiltonian H_ ') it follows 
that 

H^\AiB n ^ = E^[AlB n ^] 

(n) ~, (n) 

It is clear that the normalized cigenfunctions of the Hamiltonians H_ and H_ are related by 

^ k n \x) = (E k - Ej-^BlAnMx)] and ^> k {x) - (E k - J B n )- 1 [Aj l B„^™ ) (x)] 

with k ^ n. The operators A n or B n destroy a node in the eigenfunctions, but they are followed respectively by the 
operators B\ or A^ n that create an extra node. So the overall number of the nodes does not change. In addition, 



the normalization docs not require positive semi-definitcness of the energy eigenvalues, as in the standard case. The 
operator B^A n can be factored and compared with the above mentioned operator A = d 2 + rj{x)d x + 7(2;). One 

realizes that — r]{x) = f n (x) and —7(1') = V_ (x) + f n {x)W n (x). 

For any n there is always one missing state k = n (since A n ip n {x) = 0) which can be obtained by solving the first 

~(n) 

order differential equation B n ^\) n = 0: 



d^\x) 
dx 



7- w+ ^ws «-'"-s 



In 



t,. 



C + f x Jl{s)ds 



^\x) 



In other words, 



fi»\x) = N{C)x 



t,, 



C + f*J 2 (s)ds 



£>(*») 



£(") 



with the corresponding energy E„ = 0. All other energy eigenvalues satisfy E k = Ek — E„. Negative energy 
states appear when n > 0. The normalization constant N(C) depends on the parameter C and other parameters of 
the initial potential V_ (x). It is a constraint that allows one to determine the values of C for which the potentials 



V 



(n) 



(x) are nonsingular and cigenfunctions ipj™ ( x ) ar e well-defined. 



Let us apply the developed theory to the harmonic oscillator potential (with ui = 2) whose wave functions tjjk(x) = 
Hk(x)e~ x I" 1 and energy eigenvalues Ek — 2k arc well known. Considering first the n = case, one obtains the 
following partner potentials 



Vl°\x) =x 2 -l and V m (x) = x 2 



'dx \C + J*e~ s2 ds 



The latter is free of singularities when \C\ > y/n/2 (as follows from the normalization constant N(C)) and asymptot- 
ically behaves like x 2 — 1 (see Fig. 1). 





FIG. 1: Plot of potential V_ (x) with C = 1 and unnor- 
malized probability densities (dashed line at the corre- 
sponding level position) for its three lowest energy eigen- 
values. The limit C — > 00 corresponds to the potential 
Vi 0) (a;) = x 2 - 1 (thick blue line). 



FIG. 2: Plot of potential V_{x) with C = 3 and unnor- 
malized probability densities (dashed line at the corre- 
sponding level position) for its three lowest energy eigen- 
values. The limit C — > 00 corresponds to the potential 
V^\x) = x 2 - 3 (thick blue line). 



Its eigenvalues are E k = 2k, that is identical to those of the harmonic oscillator and the normalized ground state 



wave function reads 



Uc 2 - 7T e~ x I 2 
4 ^ C + J^e- 



' ds 



All other wave functions are derived through BqAqiP^ or by applying the fifth order raising operator (BqAq)Aq(AqBq) 
repeatedly to the normalized first excited state wave function of the potential V_ (x) : 

r(0)/^ _ ( 4 V /4 er Zx2 l 2 {\ + 2Cxe x °~ + 0Fxerf(x)e x2 



cw 



7r/ 2C + v / 7rerf(x) 



The same result was obtained by Mielnik [9| and using the Gelfand-Levitan formalism by Abraham and Moses [T(j] ■ 
Mielnik pointed out the third order raising operator for this problem b'a'b where b and tf convert the eigenf unctions 
of the partner Hamiltonians into one another and a* 1 is the usual creation operator. We immediately recognize that 
b = AqBq and a^ = A - this is just the consequence of non-uniqueness of factorization (and in the first instance of 
double supersymmetry) . Being able to reproduce this result is noteworthy, but another important fact is that one can 
obtain infinitely many such potentials (since the harmonic oscillator has infinitely many wave functions) and Miclnik's 
result is just one of them, namely n — 0. We show one more example and then turn to the most general case. This 
time letting n — 1 we find a pair of partner potentials: 



„2 

V^ (x) = x 2 - 3 and V^> (x) = x 



(l)^_„2 _o a „j f/(l)^_ Jl ~ -«' Xe 




C-2xe~ x2 +v^Ferf(x) 



The potential V_ (x) is nonsingular for \C\ > y/n (as follows from normalizing ip[ (x)) and behaves like x 2 — 3 for 

large x (see Fig. 2). Its eigenvalues are E k = 2(k — 1), i.e. the same as for V_ (x) and there is one negative energy 
state as expected since we used the first excited state of the original potential. The normalized first excited state oi 
the potential V_ (x) corresponds to the energy E{ ' — and satisfies B\^)\ (x) = 0: 



7(i), , IC 2 -7T 2xe~ x I 2 

m w = 



2^ C-2xe- :r2 + 0Ferf(x)' 
The normalized lowest state wave function is given by 

.IM, X) riV /4 e- 2 / 2 (C + 0Ferf(x)) 
Vo [ ' \ir) Ce* 2 -2x + V^erf(x)e* 2 ' 

Other states k = 2, 3, • • • are obtained by the usual construction or by applying the raising operator (B\Ai)Aq(A\B\) 

to the wave function ^ ( x )- 

At least in the case of the harmonic oscillator one can provide explicit expressions for the pairs of isospectral 
potentials generated by an arbitrary excited state ip n (x). To this end note that we need to evaluate the integral (with 
xq = 0) that appears in the denominator of the function f n {x): 



71 om 

e- 8 -Hl{s)ds = dse-*- £ — 



o Jo Jo 



m=0 



(n — m)\ 



H2(n-m)(s) 



^ (n-m)\ (to 3 V 2 2 



m=0 



-x 2 



where iFi{a; j3; z) is the confluent hypergeometric function (cf. eq. (46) in ref. [7|). 

Up until now we were assuming that for certain values of the parameter C the function f n (x) is free of singularities. 
Of course, / n (x) is non-singular if its denominator is nodclcss and latter increases monotonically (the derivative of 
the denominator is positive for all x). So the requirement is satisfied if the limits 



lim <C+ ^l(s)ds 

x^ioo I J Q 



both have the same sign, i.e. when 

|C|> 



o 

and this is true for normalizable wave functions of any potential. In our case we find \C\ > 2 n_1 n!y / 7r which is in 
agreement with the above considered particular cases n = and n = 1. Thus, one obtains non-singular potentials 
V_(x) strictly isospectral to the harmonic oscillator V_ = x 2 — (2n + 1). The same result (up to the choice of the 
energy scale) was explored for the first time in ref. |7j . We would like to stress that for a given n we supersymmetrize 
the original potential by subtracting the appropriate energy. 

Continuing our theme of potentials with the oscillator spectrum let us finally apply the presented method to 
the exactly solvable non-polynomial one-dimensional quantum potential discovered recently by Carinena, Pcrelomov, 
Rahada and Santandcr (CPRS) [llj. In their original manuscript the authors solved the Schrrodinger equation for 
the potential 

V C rps{x) = tt+4- 



2 (2:r 2 + l) 2 

assuming that the solution can be expressed as a power series in x and deriving the recursion formula which determines 
the coefficients of the expansion. Later a much more simple solution was obtained by Fellows and Smith [12J who 

showed that the potential V_ (x) = 2Vcprs + 3 is a supersymmetric partner of the harmonic oscillator x 2 + 5. We 
further analyze the CPRS potential with the aim to find new potentials with the spectrum of a simple oscillator. The 
energy eigenvalues and unnormalized wave functions of the potential V_ (x) are given by E k = 2fc and 



H k {x) + AkH k _ 2 {x) + 4fc(fc - 3)ff fc - 4 (30 „ x 2 /2 = P k (x) p _ x * /2 
2x 2 + 1 6 2a; 2 + 1 " 



M*) = ^ T1 e / 

with k = 0, 3, 4, 5, . . . Knowing one of the excited state wave functions, say 

2x z + 1 
one can generate a pair of potentials with identical spectra (see Fig. 3): 



1 |C = 21.5 



FIG. 3: Plot of the potential V_ (x) with C = 21.5 (close to 12^/tF) and unnormalized probability densities (dashed line at 
the corresponding level position) for its three lowest energy eigenvalues. The limit C — > oo corresponds to the CPRS potential 
Vi 3) (:r) (thick blue line). 



The latter is free of singularities for \C\ > 12^/tt. It might seem very unexpected that this highly non-trivial potential 
possesses such a simple spectrum as E k = 2k — 6 for k = 0, 3, 4, 5, . . . Its first excited state wave function reads 



TO), x C 2 -U4ir Ax(2x 2 + 3)e x / 2 

^3 0) = 



240? e K2 (2x 2 + l)(C + 12V^erf(a;)) - 8x(2x 2 + 3) ' 



The eigenfunctions with k ^ n are obtained by applying the operator B^A n to the wave functions ipk{x) of the original 

potential V_ (x). 

It turns out that for this problem too one can write the partner potentials generated by an arbitrary excited state 
tp n {x). We need the following two identities which are valid for any n > 3 (proofs can be found in [111]): 



P„(x)e x d t i? n _ 3 (x)e" 



J* (a:) = 4n(2a; 2 + l)H n _ 3 (x) and "V, T^ = ~ 2 ^T 



(2x 2 + l) 2 dx \ 2a; 2 + 1 J ' 

Using these identities and integrating by parts we find 

a( W 2P n {x)H n ^{x)er* 2 r _, a 2 

tp n (s)ds = +8n/ e H n _ 3 {s)ds 

tx -+- 1 J 

2P»(x)g»_ 3 (x)e- 2 2 ^ (-1)" 2"— 3 (2m)! / 13 

—5 hx Bn((n — 3)!) > — —7-7 — p^— ii<i m + - - — x 

2x 2 + l u 1 > £^1 n-m-3)! (to!) 3 V 2 2 

m— v 



Hence for any given n, the potentials 

v - l){x) = W)= x2 + & 0TW ~ 2n + 3 and V{U){x) = V - n){x) ~ V ' n{x) 

are strictly isospectral with energy eigenvalues 2(fc — n) where k, n = 0, 3, 4, 5, The potentials V_ (a;) arc non- 
singular for I CI > (n — 3)!2™ _1 ny / 7r with n > 3. We calculated V_ (a;) using the following identity which does not 
appear in ref. [ll| where some properties of the polynomials P n (x) have been studied: 

2x(2x 2 + 5)P; t (aO _ P%W = ^ 
2a; 2 + 1 P n (x) P n (x) 

The canonical form for this identity is the Schrodinger equation for the CPRS potential. 

The above discussion illustrates a general method for constructing isospectral potentials in the second-order SUSY 
QM. Given the wave functions ipu and energy eigenvalues E}. — Eq of the Hamiltonian H_ we automatically get 
the wave functions ipk and energy eigenvalues Ek — E n of the Hamiltonians H _ . This allows one to construct 
the Hamiltonian H_ whose energy spectra is identical to that of H_ and eigenfunctions arc given by rp^ 1 ' (a;) ex 
B^Anipkix) for k =^ n. Intertwining operators of higher order in derivatives have been studied [1J|. For the third-order 
reducible operator there will be two atypical Hamiltonians arising from the non-uniqueness of factorization. If the 
possibility of using operators is available (as when one uses excited states) , then the operators inherit the property of 
being ladder operators at each stage. 

We presented simple examples to show the technique of construc ting isospectral potentials, and we have tested our 



results on the consistency with other methods. It is well-known 14| that ladder operators can be constructed for 
other potentials such as Morse, Poschl- Teller, and the infinite square well. The raising and lowering operators for 
the doublc-supcrsymmctric partner potentials can be obtained along the same lines as illustrated in this paper. This 
avoids the unnecessarily technical challenges in the general problem, and these potentials will be the subject of future 
investigations. 
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